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DNA  (ssDNA)  from  a  complex  mixture  of  DNA,  containing  strands  of  different  composition  and  size.  We  have  achieved  both 
goals  and,  in  some  cases,  expanded  the  investigation  from  the  original  tasks. 

The  results  of  our  studies,  elaborated  in  the  following  sections,  suggest  that  lab-on-chip  nanofluidic  platforms  may  enable  rapid 
and  inexpensive,  characterization  and  analysis  of  DNA  biomarkers.  Advantages  include  overall  ease  of  operation  of  the  device, 
which  does  not  require  loading  of  gels  as  sieving  matrices;  inexpensive  analysis,  as  the  only  reagents  required  are  water-based 
buffers;  fast  identification,  since  sample  signatures  can  be  produced  in  the  order  of  less  than  a  minute. 

Findings  from  this  work  have  resulted  in  peer-reviewed  publications  (one  conference  and  one  journal  paper).  Potential  target 
applications  of  our  identification  platform  include  analysis  of  DNA  biomarkers  for  environmental  pollutants  or  early  disease 
diagnosis;  analysis  of  mitochondrial  DNA  for  forensic  identification:  investigation  of  protein  kinetics.  
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Discussion  of  Final  Results 

AIM  1:  Demonstrate  FASS  for  specific  DNA  samples,  optimizing  for  buffer  concentration 
and  pH,  to  determine  the  largest  DNA  concentration  possible  in  the  shortest  time. 

Nanofluidic  technology  is  gaining  popularity  for  bioanalytical  applications  due  to  advances  in 
both  nanofabrication  and  design.  One  major  obstacle  in  the  widespread  adoption  of  such 
technology  for  bioanalytical  systems  is  efficient  detection  of  samples  due  to  the  inherently  low 
analyte  concentrations  present  in  such  systems.  This  problem  is  exacerbated  by  the  push  for 
electronic  detection,  which  requires  an  even  higher  sensor-local  sample  concentration  than 
optical  detection.  We  have  preliminary  results  showing  the  potential  of  FASS  as  an  effective  pre¬ 
concentration  mechanism  at  the  nanoscale  [1],  However,  before  this  project  we  had  not  explored 
the  case  of  samples  consisting  of  complex  mixtures  of  DNA.  Using  FASS  would  allow  for 
obtaining  an  on-board  concentration  method  that  eliminates  PCR  amplification  and  overcomes 
the  related  drawbacks  that  may  lead  to  unreliable  analysis  results,  and  prevent  absolute 
quantification  of  analytes. 

A.l  Chemicals  and  reagents 

For  all  experiments,  we  used  potassium  phosphate  buffer  at  pH~7.2,  obtained  by  monobasic  and 
dibasic  sodium  phosphate  (Sigma  Aldrich,  Inc.)  in  the  appropriate  ratios,  dissolved  in  deionized 
water  (Millipore,  Inc.  pH  7.5).  Stock  solution  of  200  mM  was  prepared  using  de-ionized  filtered 
water  at  18  MQ/cm  provided  by  a  Millipore  Milli-Q  Gradient  filtration  system.  The  stock 
solution  was  then  diluted  to  concentrations  of  0.5,  1.0,  2.5,  5.0,  10.0,  and  25  mM  and  seeded 
with  10  p  M  sodium  tluoresccin  (Sigma  Aldrich,  Inc.).  The  pH  and  conductivity  were  then 
measured  (Oakton  pH  and  conductivity  meter)  and  the  solution  was  filtered  with  0.2  pm 
polytetrafluoroethylene  syringe  filter. 

A.l  DNA  samples 

Fluorescently  labeled  (FAM  tag)  DNA  oligomers  (10,  20,  and  50  bases  long)  were  purchased 
with  standard  desalting  and  additional  HPLC  purification  for  the  50  base  strands,  from  Integrated 
DNA  Technologies,  Inc. 

A.l  Experimental  setup 

The  experimental  setup  is  diagramed  in  Fig.  1.  The  motion  and  spatial  distribution  of  the  sodium 
fluorescein  was  imaged  using  an  inverted  epifluorescent  microscope  (Olympus  1X70)  fitted  with 
a  60/  1.00  NA  (NA  denotes  numerical  aperture)  water  objective  lens  (Olympus,  Inc.). 
Illumination  was  provided  from  a  200  W  Hg-arc  lamp  filtered  both  before  and  after  the  sample 
using  interference  filters  and  a  dichroic  mirror  specific  to  the  peak  fluorescein  absorption  and 
emission  wavelengths  of  485  and  535  nm,  respectively.  A  ProScan  III  17  automated  stage  with  a 
resolution  of  0.1  m  was  used  to  monitor  the  development  of  the  sample  ions  as  they  propagate 
down  the  East  channel.  We  applied  voltages  to  each  well  of  the  nanofluidic  channel  through 
platinum  electrodes  controlled  by  a  multichannel  high  voltage  sequencer  (HVS448-6000D, 
Labsmith,  Inc.).  A  preprogrammed  sequence  of  voltages  established  the  desired  fields  for 
flushing,  sample  loading,  gating,  and  injection/separation  estimated  using  finite  element 
simulations  (Comsol,  Inc.).  We  acquired  intensity  data  using  a  back  illuminated  and  intensified 
EMCCD  camera  (Ixon+,  Andor  Co.)  with  a  512x512  pixel  array  and  16-bit  digitization.  Images 
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were  corrected  using  MATLAB  by  subtracting  a  background  image  from  the  raw  data  and 
normalizing  by  the  difference  between  the  flatfield  and  background  image  following  [23]. 


Figure  1:  Diagram  of  the  experimental  setup.  A  Prior  automated  stage  (ProScan  HI  17)  allowed  for  real¬ 
time  observation  of  the  sample  as  it  passed  through  the  detection  points  of  5,  10,  15,  20,  or  25  mm  from 
the  intersection. 


The  interface  between  high-  and  lowconductivity  buffer  regions  was  generated  with  potassium 
phosphate  buffer  solution,  driving  the  lower  conductivity  buffer  containing  the  sample  from  the 
North  to  the  South  well,  and  the  higher  conductivity  buffer  from  the  West  and  the  East  to  the 
South  well  (Fig.  2,  step  a).  Once  a  conductivity  interface  was  established,  a  gating  method  was 
necessary  to  inject  enough  low-conductivity  buffers  into  the  channel  to  allow  for  fluorescein  to 
enter  the  East  channel  (Fig.  2,  step  b).  Once  enough  fluorescein  entered  the  East  channel, 
injection  voltages  were  used  to  propagate  the  plug  down  the  channel  (Fig.  2,  step  c). 
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Figure  2:  Schematic  of  sample  injection  sequenee  with  flow  direction  denoted  with  arrows  and  the  low- 
conductivity  sample  region  is  denoted  by  the  shaded  area,  (a)  First,  the  low-conductivity  sample  is  loaded 
from  the  North  well  to  the  South  well  In  this  step,  high-conductivity  buffer  is  flowing  from  both  the  East 
and  West  wells  toward  the  South  well  (b)  Next,  gating  is  performed  for  several  seconds  in  order  to  inject 
a  long  plug  of  low-eonduetivity  containing  the  fluorcsecnt  sample.  Low-eonduetivity  fluid  and  sample 
ions  enter  the  East  channel  from  the  North  well,  (c)  Finally,  the  sample  is  injected  by  switching  to  the 
injection  voltages.  Here,  high  conductivity  buffer  flows  from  the  West  to  the  East  channel,  closing  off  the 
introduction  of  sample  from  the  North  well  and  allowing  for  the  sample  to  propagate  as  a  plug  down  the 
channel. 


A.l  Experimental  Results 

Our  first  experiments  served  as  control  and  were  performed  with  fluorescein,  confirming  that  at 
the  nanoscale  FASS  depends  on  a  number  of  parameters,  such  as  channel  height  (Figure  3), 
electrolyte  concentration  (Figure  4),  electric  field,  gating  time,  etc. 
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Figure  3:  Comparison  of  FASS  concentration  enhancement  in  lpm  and  250nm  channels.  The  analyte  is 
lOpM  fluorescein  in  l:10mM  phosphate  buffer.  Note  that  enhancement  in  nanoehannels  can  be  much 
larger  than  in  mieroehannels  (up  to  two  orders  of  magnitude  [1]). 
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Figure  4:  Study  of  FASS  dependence  on  buffer  concentration  at  the  nanoseale  in  same  experimental 
conditions  as  in  Figure  1.  Note  that,  while  in  1  pm  channels  the  enhancement  is  roughly  constant  with 
concentration  ratio,  in  250  nra  channels  enhancement  strongly  depends  on  background  salt  concentration: 
the  lower  the  concentration  (e.g.,  larger  EDL),  the  higher  the  enhancement. 


A  qualitatively  similar  behavior  was  observed  in  experiments  with  DNA  samples  (lObp  ssDNA), 
as  shown  in  Figure  5,  in  terms  of  dependence  on  buffer  concentration,  channel  size  and  gating 
time. 


lum  FASS  lObp  ssDNA  (10:1mM  Buffer  Ratio) 
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Figure  5:  DNA  PASS  experiments  in  1  pm  (a-c)  and  250nm  (b-d-e)  channels  and  different  buffer 
conditions:  (a-b)  10:lmM,  (c-d)  20:lmM;  (e)  different  gating  times  in  10:lmM  buffer.  DNA  is  lObp  ss, 
as  described  in  materials  section.  Error  bars  represent  one  standard  deviation. 

We  have  observed  that,  for  DNA  injections,  the  gating  time  was  a  more  critical  parameter  than 
fluorescein,  in  terms  of  obtaining  consistent  sample  plugs.  The  inconsistency  results  in  a 
different  amount  of  sample  entering  the  injection  channel,  with  an  obvious  impact  on  the 
maximum  obtainable  concentration  enhancement.  In  order  to  better  understand  the  origin  of  this 
phenomenon  and  optimize  injections,  we  have  developed  a  COMSOL  model,  consisting  of  three 
systems  of  coupled  equations:  1)  the  mass  transport  equation, 

dCj 

—  +  V  •  (— DVq  -  ZiitjniFctV V)  +  u  •  Vct-  =  0 


describing  the  transport  of  (two)  dilute  species,  c,-,  with  electromigration  and  convection  terms:  cy 
which  models  the  buffer  concentration  throughout  the  channel,  and  C2  which  represents  a  charged 
species  (e.g.,  DNA). 


2)  Gauss  law,  which  describes  the  electric  currents 

-V  ■  (eW)  =  0 

and,  thus,  governs  the  electric  potential  in  a  conductive  media,  which  we  couple  to  the  buffer 
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concentration  to  simulate  conductivity  differences  between  buffers  of  different  concentration. 
3)  the  continuity  and  Navier-Stokes  equations 

V  •  u  =  0 

V  ■  [-pi  +  p  (Vu  +  (Vu)T)]  =  0 

which  govern  fluid  flow  within  the  channel  with  electroosmotic  flow  boundary  conditions: 


U  —  flEO^t 


which  creates  wall  velocity  boundary  conditions  equal  to  the  electroomotic  mobility  of  the 
channel.  This  is  coupled  to  the  local  electric  field  from  the  electric  currents  module,  as  well  as  a 
zeta  potential,  which  is  dependent  on  the  local  buffer  concentration. 

Figure  6  shows  a  parametric  study  of  the  time  it  takes  for  fluorescein  to  enter  the  injection 
channel  for  different  buffer  concentration  ratios.  We  note  that,  given  an  electric  field,  for 
increasing  buffer  ratios,  the  time  fluorescein  takes  to  “enter”  the  channel  (reach  a  location  few 
microns  from  the  intersection)  during  gating  is  longer.  During  this  time,  unstable  fluid  behavior 
at  the  intersection  prevents  fluorescein  from  entering.  The  dashed  line,  t2,  describes  the  time  for 
the  transient  unstable  behavior  to  disappear,  and  a  steady  state  behavior  to  emerge. 


Electric  Field  [kV/m] 


Figure  6:  Experimental  study  of  instabilities  in  FASS  for  different  buffer  concentration  ratios:  l:lmM 
(blue  -  control),  2:lmM  (green),  5:lmM  (red),  10:lmM  (black),  and  50:lmM  (purple).  In  all  cases,  the 
solid  line  represents  the  time  ti  when  the  sample  first  enters  the  east  (injection)  channel,  while  the  dashed 
line  represents  the  time  t2  when  sample  reaches  the  end  of  the  field  of  view  on  the  east  channel. 

In  particular,  at  lower  buffer  ratios,  the  sample  is  let  in  right  away  (Figure  7  a-b).  At  higher 
buffer  ratios,  sample  is  “concentrated”  right  at  the  junction  before  actually  entering  the  injection 
channel  (Figure  7  d).  The  5:1  ratio  seems  to  be  the  “middle”  case  (Figure  7  c):  after  gating, 
sample  is  let  into  the  injection  channel  (data  not  shown),  then  it  is  retracted  back  at  the  junction 
where  it  focuses.  The  same  figure  shows  the  very  good  qualitative  agreement  obtained  by  our 
numerical  simulation. 
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10:1  GVo  E-field  ;  Time  =  55ms 


Figure  7:  False  color  images  of  injection  experiments  (a-d)  and  numerical  simulations  (e-h)  for  different 
buffer  ratios,  and  same  electric  field  (GVo  -  120kV/m). 


Gating  Time  [sec] 


Figure  8:  Experimental  results  showing  plug  length  as  a  function  of  buffer  concentration  for  different 
electric  fields:  1)  120kV/m,  b)  60kV/m,  and  c)  15kV/m. 


Figure  8  shows  experimental  data  relating  plug  length  to  electric  field  strength,  for  different 
buffer  concentration  ratios.  As  expected,  for  lower  electric  fields  and  gating  times,  the  sample 
plugs  are  shorter.  However,  we  notice  the  nonlinear  behavior  for  higher  buffer  concentration 
ratios,  which  indicates  unstable  fluid  behavior  at  the  intersection.  These  flow  instabilities  prevent 
fluorescein  from  entering  the  channel,  resulting  in  lower  amount  of  sample  injected  and  reduced 
concentration  enhancement.  These  and  other  experimental  data  will  be  used  to  improve  the 
model  so  that  it  can  be  used  to  predict  future  experiments. 
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Figure  9:  Concentration  enhancement  for  DNA  (red)  and  fluorecein  (blue)  in:  a)  lpm  and  b)  250nm 
channels  for  different  buffer  concentration  ratios. 

Figure  9  presents  PASS  experiments  results  for  fluorescein  (blue)  and  DNA  (red)  in  1pm  and 
250nm  channels,  for  different  buffer  concentration  ratios.  For  DNA  and  fluorescein  in  the 
microchannel  case  (Fig.  9a),  the  10:1  and  5:0. 5mM  cases  yield  approximately  the  same 
enhancement,  equal  to  the  conductivity  ratio,  although  the  DNA  values  are  consistently  higher. 
However,  for  the  20:lmM  case,  DNA  FASS  results  in  a  much  higher  enhancement,  which  we 
postulate  might  be  due  to  the  higher  negative  charge  of  DNA.  In  fact,  from  [  1  ]  we  know  that  also 
in  microchannels  the  maximum  enhancement  non-linearly  depends  not  only  on  the  EDL 
thickness,  but  also  on  the  analyte  charge  (see  Fig.  10  [1],  reported  below  for  completeness). 
Moreover,  the  transition  to  the  exponential  growth  region  in  enhancement  strongly  depends  on 
the  charge  of  the  analyte.  Taking  the  DNA  electrophoretic  mobility  pep  =  -3.9e"*  [m2/Vs]  and 
diffusivity  D  =  2e"10  [m2/s],  we  can  roughly  estimate  the  charge  number  of  a  lOnt  DNA  to  be  z  = 
~3:  preliminary  calculations,  using  [1],  predict  a  behavior  for  consistent  with  measured  data  (Fig. 
9a). 


Figure  10  [1]:  Maximum  theoretical  concentration  enhancement  ratio  as  a  function  of  Debye  length  for  a 
5  mm  long  plug  in  a  30  mm  straight  channel,  plotted  along  with  experimental  data.  In  this  figure  the 
Debye  length  is  defined  for  the  low-conductivity  region  and  the  conductivity  ratio  is  10  to  1.  For  this 
calculation,  the  concentrations  arc  held  fixed  at  1  and  10  mM  and  the  channel  height  is  varied.  The  250 
nm  data  are  shown  as  circles  and  the  1  pm  channel  data  are  shown  as  triangles.  The  experimental  data  are 
below  the  theoretical  maximum  due  to  diffusion  and  other  inaccuracies  associated  with  gating  and  plug 
length,  however  the  trend  is  in  good  agreement  with  the  data.  The  open  shapes  are  the  lower  electric  field 
(16  kV/m)  and  the  filled  points  are  for  the  higher  electric  fields  (32  kV/m). 
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Regarding  the  250nm  channels  (Fig.  9b),  fluorescein  data  shows  the  expected  behavior, 
consistent  with  past  results,  concentration  enhancement  is  much  lower  than  expected.  To 
understand  enhancement  limits  in  250nm  channels,  we  started  by  performing  experiments  with 
different  gating  times:  since  the  concentration  enhancement  depends  on  amount  injected  (which 
will  be  different  at  different  buffer  ratios),  are  we  injecting  enough  DNA?  Or  are  electrostatic 
interactions  between  DNA  strands  resulting  in  less  analyte  being  injected? 

At  first,  we  recorded  and  measured  the  shape  and  area  under  the  electropherograms  at  different 
observation  points,  to  see  the  progress  of  analyte  stacking.  To  our  surprise  (Figure  1 1),  we  found 
that,  instead  of  remaining  constant  (conservation  of  mass),  the  area  under  the  electropherograms 
was  increasing  at  later  observation  points.  Possible  explanations  for  such  behavior  are:  1)  the 
quantum  efficiency  of  the  FAM  dye,  which  is  pH  dependent,  might  cause  an  increase  in 
brightness  as  the  DNA  enters  the  high  conductivity  region;  2)  while  traveling  towards  the  East 
channel,  the  sample  plug  slows  down,  thus  resulting  in  peak  broadening  and  an  apparent  larger 
area. 

A  new  sample  batch  consisting  of  DNA  tagged  with  Alexa  448  has  been  ordered  and  will  be 
used  to  continue  this  investigation  after  the  end  of  this  project.  Preliminary  measurements  of 
plug  speed  within  250nm  channels  do  not  show  any  appreciable  difference  (i.e.,  slowing  down) 
of  DNA  compared  to  fluorescein  and  compared  to  what  was  expected  based  on  our  model  [1], 
More  detailed  studies  are  needed  to  understand  the  origin  of  this  observed  phenomenon. 


lum  lObp  saDNA  (20  ImM  Buffer  Ratio)  16kV/m 


Time  Is] 


Figure  11:  Example  of  FASS  DNA  (lOnt  ssDNA,  20:lmM  buffer,  16kV  electric  field  injection) 
electropherograms  collected  at  different  observation  points.  Notice  how  the  area  under  each  curve,  instead 
of  remaining  approximately  constant,  increases  with  time. 

A.l  Conclusions 

We  have  demonstrated  that  FASS  is  a  viable  technique  to  obtain  ob-board  concentration 
enhancement  of  DNA,  in  both  micro  and  nanochannels.  As  part  of  this  work,  we  have  also 
developed  a  COMSOL  prototype  model  that  will  allow  for  studying  the  instability  behavior 
during  injection:  such  model  will  prove  useful  in  tuning  injection  parameters  (such  as  electric 
fields  as  a  function  of  buffer  concentrations)  to  obtain  repeatable  and  efficient  injections. 
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Our  experiments  showed  that  DNA  does  not  behave  like  fluorescein  and  that,  therefore,  the 
Matlab  model  we  had  previously  developed  needs  to  be  modified.  The  origin  of  the  discrepancy 
in  behavior  are  not  yet  clear,  and  further  experiments  are  needed. 

AIM  2:  Demonstrate  identification  of  target(s)  ssDNA  from  a  complex  mixture  of  DNA, 
containing  strands  of  different  composition  and  size. 

We  started  our  investigation  by  proving  that,  in  nanofluidic  channels,  different  ssDNA  could  be 
separated  from  each  other  and  from  the  respective  dsDNA  obtained  by  hybridization  with  the 
corresponding  complement  (Figure  12  a-b  and  Figure  14).  Figure  12  e)  also  shows  that  nanoscale 
injections  can  yield  information  on  sample  concentration,  which  is  proportional  to  the  area  under 
the  peak  on  the  corresponding  eleetropherogram.  However,  in  analyzing  the  data,  we  found  that 
in  addition  to  the  expected  changes  in  mobility  with  buffer  ionic  strength  and  channel  size,  the 
relative  changes  in  mobility  between  ss  and  dsDNA  become  more  pronounced  as  ehannel  size 
decreases  (see  Figure  12  e),  and  that  while  confinement  seems  to  dominate  mobility  for  ssDNA, 
dsDNA  seem  to  be  more  affected  by  changes  in  EDL.  These  results  posed  fundamental  questions 
on  the  dynamics  of  short  DNA  strands  in  micro  and  nano-eonfinements,  whieh  a  survey  of 
published  literature  could  not  answer.  In  fact,  although  there  has  been  significant  progress  in  the 
experimental  characterization  of  micro  and  nanotluidic  bioanalytical  platforms  [5-8],  the  precise 
mechanism  involved  with  separations  of  biomolecules  at  the  micro  and  nanoscale  is  not  fully 
understood.  This  is  due  in  part  to  the  nature  of  biomolecules,  which  can  have  complex 
interactions  with  each  other,  as  well  as  with  their  environment. 

To  investigate  the  behavior  of  short-stranded  DNA  in  micro  and  nano  confinements,  we  have 
performed  a  parametric  experimental  study  on  the  transport  and  separation  behavior  of  ss-  and 
dsDNA  (10,  20,  and  50  bases  long),  under  electrokinetic  flows  in  micro-  and  nanochannels.  We 
hypothesize  that  several  faetors  impact  electrophoretic  mobility:  (i)  hybridization  interactions, 
(ii)  salt-dependent  electric  double  layer  (EDL)  changes  around  the  DNA,  which  alter  the 
size/conformation,  (iii)  differences  in  drag  due  to  length  dependent  conformational  differences, 
(iv)  hydrodynamic  confinement  effects  due  to  decreased  channel  size,  and  (v)  electrostatic 
repulsion  from  the  walls.  These  factors  have  been  studied  using  both  analytical  and  experimental 
observations  towards  the  optimization  of  an  effective  platform  for  short  DNA  separation  and 
identification. 

A.2  Chemicals  and  reagents: 

Buffered  solutions  consisted  of  a  3:17  ratio  of  monobasic  and  dibasie  sodium  phosphate, 
dissolved  in  deionized  water  (Millipore,  Inc.  pH  7.5)  at  10  mM  concentration  (1.5  mM 
monosodium  phosphate,  8.5  mM  disodium  phosphate).  Buffer  ionic  strength  was  varied  by  the 
addition  of  NaCl  in  concentrations  from  5  to  50  mM.  All  solutions  were  filtered  with  200  nm 
pore  syringe  filters  (Nalgene,  Rochester,  NY)  prior  to  use. 

A.2  DNA  samples: 

DNA  oligomers  (10,  20,  50  nt  long)  were  purchased  with  standard  desalting  and  additional 
HPLC  purification  for  the  50  base  strands,  from  Integrated  DNA  Technologies,  Inc.  Sequences 
were  designed  to  be  a  repetition  of  the  same  basic  1 0  nucleotide  “unit”  (5'-AA  GAG  GAG  GG  - 
3'),  so  that  changes  in  electrophoretic  mobility  eould  be  attributed  to  differences  in  length 
dependent  conformations  rather  than  differences  due  to  sequence  or  self-complementarity.  For 
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each  sequence,  we  use  both  a  fluorescently  labeled,  5’  modified  (6-FAM)  strand  and  its  untagged 
complement. 


10  20  30  40  50  60 

Ionic  Strength  [mM] 


lOOnm  ss 
2  50  run  ss 
lum  ss 


x10" 

- r 

-3 


b) 


-fr-IOOnm  ds 

- 250nm  ds 

-*-1  um  ds 


4^0  20  30  40  50  60 

Ionic  Strength  [mM] 


Figure  12:  a)  and  b)  Electrophoretic  mobility  measurements  of,  respectively,  10  ss-  and  dsDNA 
measured  in  a  1  pm  (blue  trace),  250nm  (red  trace)  and  100  nm  (green  trace)  deep  channels  at  different 
ionie  strengths.  Error  bars  correspond  to  the  standard  deviation  of  all  trials  for  a  given  experiment;  c) 
same  data  but  normalized,  for  each  channel  depth,  by  subtracting  the  corresponding  mobility  value  at 
lOmM  ionic  strength;  d)  same  data  but  normalized  by  dividing  the  mobility  values  at  each  ionic  strength 
by  the  corresponding  value  for  1  pm  deep  ehanncl;  e)  quantification  of  DNA  concentration  (50  nt  ssDNA 
in  lOmM  phosphate  buffer)  in  a  nanochannel. 

A.2  Microfluidic  device  and  setup: 

MicroChannel  devices  were  designed  and  custom  fabricated  in  borosilicate  wafers  using 
conventional  MEMS  processing  techniques  (Dolomite  Ltd,  UK).  The  devices  consist  of  a  simple 
cross  channel  geometry,  either  100  nm  deep  and  7  pm  wide,  or  1  pm  deep  and  9  pm  wide,  with 
four  inlet/outlet  ports  (labeled  N,  S,  E,  W,  shown  in  Figure  13).  Electrical  potentials  were  applied 


10 


Pennathur,  UCSB  W9132T-12-2-0023 


at  these  wells  using  platinum  electrodes  (Omega  Eng.  Inc.,  Stamford,  CT),  connected  to  a  high 
voltage  power  supply  (LabSmith  HVS448).  The  pre-programmed  voltage  scheme  for  sample 
loading  and  injection  was  designed  following  the  work  of  [19].  In  the  loading  step,  voltages  are 
set  such  that  the  fluid  flows  from  top  to  bottom  (N-S)  in  the  vertical  channels  and  toward  the 
intersection  in  the  horizontal  channels.  In  the  sample  injection  step,  the  applied  voltages  are 
instantaneously  switched  (W  at  the  highest  voltage,  N,  S  at  high,  E  at  ground),  and  the  sample  is 
injected  in  the  horizontal  channel  under  an  applied  electric  field  of  16kV/m.  During  the  sample 
injection  step,  fluid  flows  from  left  to  right  in  the  horizontal  channel  and  away  from  the 
intersection  in  the  vertical  channels  [20].  Between  injections,  the  channel  is  flushed  using  the 
same  scheme  as  the  loading  step. 

A. 2  Fluorescence  Imaging: 

The  transport  and  spatial  distribution  of  fluorescent  analytes  were  imaged  with  an  inverted 
epifluorescence  microscope  (Olympus  1X70,  Olympus,  Inc.)  fitted  with  a  60X  water  immersion 
objective  lens  (1.0  NA,  Olympus,  Inc.)  and  a  0.5X  lens  tube,  to  increase  sensitivity.  Illumination 
from  a  200  W  Hg-arc  lamp  was  filtered  with  a  FITC  fluorescence  filter  cube  (Omega,  Inc.) 
containing  excitation  and  emission  filters  and  a  dichroic  mirror  matched  to  the  fluorescence 
spectrum  of  6-FAM.  Images  were  recorded  using  a  back  illuminated  EMCCD  camera  (Ixon 
Ultra897,  Andor  Technology)  with  a512x512  pixel  array  and  16-bit  digitization  to  a  PC  [21]. 
Frame  rate  (10  to  20  Hz)  and  exposure  time  (0.05  to  0.1  s)  varied  depending  on  the  channel 
depth  and  analyte.  Background  subtraction  was  performed  on  all  images  using  custom  Matlab 
(The  MathWorks)  programs,  to  enhance  signal  to  noise  ratio  [21]. 

A.2  Microfluidic  Electrophoretic  Mobility  Measurements: 

The  electrophoretic  mobility  was  determined  directly  from  elution  time  measurements  during 
electrokinetic  injections,  as  described  by  Pennathur,  et  al.  [20].  Briefly,  the  microscope  objective 
was  positioned  downstream  of  the  injection  point,  and  the  velocity,  vmcas,  was  computed  as  the 
ratio  of  distance  from  injection,  L,  and  the  measured  time,  t,  required  for  the  analyte  to  reach  the 
center  of  the  field  of  view  (see  Figure  13).  The  C,  potential  needed  to  compute  vEOf  from  the 
Helmholtz-Smoluchowski  equation  [22]  was  found  for  each  experiment  from  the  velocity 
measurement  of  a  fluorescein  marker  that  was  injected  with  the  DNA  sample.  The 
electrophoretic  mobility  of  fluorescein  was  determined  in  a  separate  set  of  experiments  [23], 
Although  the  calculation  of  the  electroosmotic  velocity  vE0F  is  trivial  for  microchannels,  to  find 
the  correct  electroosmotic  flow  in  nanochannels  we  numerically  solved  the  non-linear  Poisson- 
Boltzmann  equation,  and  used  the  mean  of  the  resultant  velocity  profile  [24],  Once  Veof  was 
determined,  it  was  subtracted  from  the  measured  velocity  to  give  the  electrophoretic  velocity  of 
the  analyte,  vEP.  Finally,  to  determine  electrophoretic  mobility,  pep,  the  electrophoretic  velocity 
vEP  was  divided  by  the  electric  field  E: 

H-EP  =  £  VEP  =  £  ( vmeas  ~  vEOF )  =  £  (“  —  vEOF )  (1) 
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Figure  13:  Schematic  of  sample  injection  and  separation  within  a  nanofluidic  platform,  as  well  as  an 
example  of  a  typical  electropherogram  obtained  in  our  study.  After  loading  (a),  the  sample  is  injected 
along  the  horizontal  channel  (b).  A  detector  senses  the  arrival  of  the  sample  as  a  variation  in  the  measured 
background  signal,  generating  an  electropherogram.  The  characteristic  electropherogram  shown  here  (c) 
demonstrates  the  separation  of  10  bp  ss-  (peak  2)  and  dsDNA  (peak  3)  in  a  1  pm  deep  channel; 
fluorescein  (peak  1),  is  a  mobility  marker. 

A.2  Experimental  Results  And  Discussion 

Figure  14  shows  the  electropherograms  of  10,  20,  50  bp  ss-  and  dsDNA  in  a  1  pm  channel  at 
buffer  ionic  strengths  of  10,  15,  20,  30  and  60  mM.  In  each  plot,  the  time  axis  is  scaled  to  ease 
comparisons,  according  to  the  following  formula: 

j.  _  j.  £^^1077 sr\  \ 

'-scaled, XmM  ~  '-meas,XmM  f  (AJ 

CFL.XmM 

where  tmeas,xmM  is  the  time  measured  at  ionic  strength  XmM  (i.e.,  current  experiment),  tFL,iomM 
and  tFL,xmM  are  the  fluorescein  arrival  times  at,  respectively,  lOmM  and  XmM.  Sealing 
compensates  for  the  longer  detection  times  at  higher  ionic  strengths,  which  have  lower  channel  L, 
potentials  and  therefore  slower  electroosmotic  flow  [29],  thus  allowing  for  comparison  across  all 
experimental  conditions.  We  note  that  some  electropherograms  show  more  than  the  expected 
three  peaks  (i.e.,  fluorescein,  ssDNA,  dsDNA):  we  attribute  them  to  populations  of  fragmented 
DNA,  resulting  from  the  DNA  synthesis,  and  often  observed  in  DNA  separation  experiments 
under  similar  conditions  [5],  [30]. 

DNA  behavior  and  separation  dynamics  are  affected  by  a  number  of  coupled  physical 
mechanisms,  including  (i)  hydrodynamic  interactions,  (ii)  salt-dependent  EDL  changes  around 
the  DNA,  which  alter  its  size/conformation,  (iii)  differences  in  drag  due  to  length  dependent 
conformational  differences,  (iv)  confinement  effects  due  to  shrinking  channel  size,  and  (v) 
electrostatic  repulsion  from  walls.  The  result  is  a  highly  complex  system,  where  each  parameter 
needs  to  be  carefully  designed  in  order  to  maximize  separation  resolution.  Figure  14  shows  some 
interesting  general  trends.  First,  we  note  that,  for  any  strand  length,  the  peak  separation 
resolution  increases  with  ionic  strength.  For  longer  strands,  this  effect  can  be  attributed  in  part  to 
the  longer  arrival  time  to  the  detector.  In  addition,  at  higher  ionic  strengths,  the  50  base  ssDNA 
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conformation  will  be  more  compact,  thereby  reducing  the  drag  and  lowering  the  absolute 
electrophoretic  mobility  (i.e.,  allowing  for  the  DNA  to  elute  faster). 


50  Bases 


Figure  14:  Scaled  elcctrophcrograms  of  all  ss-  and  dsDNA  experiments  in  a  1  pm  channel. 
Electropherograms  are  scaled  due  to  the  large  ionic  strength  based  changes  in  electroosmotic  mobility. 
The  scaling  aligns  the  fluorescein  peaks  to  show  the  change  in  resolution  as  a  function  of  ionic  strength. 
Melting  temperatures,  Tm,  for  DNA  at  each  condition  are  superposed  on  the  electropherograms.  Each 
electropherogram  is  plotted  as  an  average  of  raw  data  trials. 


Figure  14  also  shows  low  separation  resolution  of  the  10  bp  strands  at  low  ionic  strength.  This  is 
a  consequence  of  hybridization  interactions  (effect  (i)),  because  and  the  melting  temperature,  Tm 
(and  therefore  Kd,  the  dissociation  constant  that  describes  the  equilibrium  ratio  of  ds-  to  ssDNA), 
is  close  to  room  temperature.  As  a  result,  during  the  experiments,  the  short  complementary  DNA 
are  in  a  dynamic  equilibrium,  continuously  associating  and  disassociating,  with  significant 
impact  on  the  separation  resolution.  To  better  understand  this  dynamic,  we  computed  equilibrium 
ratios  using  the  DINAmelt  Webserver,  given  DNA  concentrations,  sequence,  and  ionic  strengths. 
From  these  values,  we  determined  initial  Kd’s.  Note  that,  due  to  the  non-equilibrium  nature  of 
our  injection  experiments,  equilibrium  ratios  are  constantly  changing  during  experiments.  For 
example,  the  equilibrium  ratio  of  ss-  to  dsDNA  for  10  base  in  lOmM  buffer  is  initially  0.0474, 
when  there  is  an  excess  of  primary  ssDNA  present.  However,  the  ratio  can  increase  to  0.1747, 
after  the  excess  ssDNA  separates  from  the  rest  of  the  plug.  Diffusion  and  dispersion  of  the 
sample  lower  the  concentration  further,  increasing  the  equilibrium  ratio  even  more.  On  Figure 
1 4,  we  write  the  range  of  Tm  calculated  from  DINAmelt  for  the  equilibrium  ratio  at  the  start  of 
the  experiment  (higher  Tm)  and  during  the  experiment  (lower  Tm).  For  a  more  detailed  discussion 
on  the  effect  of  kinetics  on  such  separations,  we  refer  the  reader  to  [3 1]. 

Finally,  we  note  the  extremely  low  signal  corresponding  to  the  50  bp  single  strand,  which 
decreases  even  further  at  higher  ionic  strengths.  A  similar  effect  can  be  seen  also  for  the  20  bp 
peak  at  the  highest  buffer  ionic  strength.  We  hypothesize  this  phenomenon  is  a  consequence  of 
using  a  FAM  tag  with  our  specific  guanine  (G)  rich  DNA  sequence,  which  has  been  shown  to 
efficiently  quench  fluorescence  [32].  Longer  ssDNA  in  high  ionic  strength  solutions  can  assume 
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more  compact  structures,  owing  to  their  flexibility  and  screening  of  electrostatic  repulsion 
between  monomers.  The  result  is  an  increase  in  the  amount  of  G  near  the  6-FAM  that  contributes 
to  quenching.  Furthermore,  since  the  50  base  strand  is  made  of  5  repeats  of  the  10  base  G-rich 
strand,  there  is  high  probability  for  G-quadruplex  formation  [33]  or  other  offset  hybridization 
configurations,  which  could  enhance  fluorescence  quenching  or  affect  the  solubility  of  the 
ssDNA.  This  effect  is  reduced  upon  hybridization  to  the  complementary  strand  to  form  stiff,  rod¬ 
like  dsDNA,  which  moves  many  of  the  G  away  from  the  6-FAM  tag  at  the  5’  end  of  the  DNA. 
Fluorescence  measurements  (not  reported  here)  of  the  6-FAM  tagged  ss-  and  dsDNA,  performed 
in  a  well  plate  using  a  Tecan  Infinite  200Pro,  confirmed  a  ~50%  decrease  in  the  fluorescence  of 
the  ssDNA  compared  to  the  dsDNA. 

The  complete  set  of  measured  values  for  the  electrophoretic  mobilities  of  ss-  and  dsDNA  in  100 
nm  and  1  pm  channels  is  shown  in  Figure  15,  as  well  as  listed  in  Table  1. 

Table  1:  Absolute  electrophoretic  mobilities  and  associated  errors  for  all  the  ss-  and  dsDNA  measured  in  1  pm  and 
100  nm  deep  microfluidie  channels.  Errors  represent  95%  confidence  intervals  according  to  the  number  of  trials 
listed  in  Table  1 . 


Absolute  Electrophoretic  Mobility  [E-8  m2/Vs] 


lOmM 

15mM 

20mM 

30  mM 

50mM 

1  pm  10  ss 

3.90  ±0.02 

3.70  ±0.01 

3.55  ±0.01 

3.37  ±  0.02 

3.07  ±  0.00 

ds 

4.12  ±0.04 

3.95  ±  0.03 

3.82  ±0.01 

3.65  ±0.01 

3.36  ±0.01 

20  ss 

3.95  ±  0.03 

3.75  ±0.03 

3.61  ±0.01 

3.44  ±0.01 

3.15  ±  0.01 

ds 

4.19  ±0.03 

3.99  ±0.02 

3.86  ±0.00 

3.69  ±0.01 

3.42  ±0.01 

50  ss 

3.96  ±0.01 

3.78  ±0.01 

3.64  ±0.01 

3.45  ±0.01 

3.15  ±  0.01 

ds 

4.24  ±0.01 

4.08  ±0.01 

3.95  ±0.00 

3.79  ±0.02 

3.38  ±0.04 

lOOnm  10  ss 

3.64  ±0.03 

3.52  ±0.02 

3.40  ±0.01 

3.24  ±0.01 

2.94 

ds 

3.81  ±0.04 

3.72  ±0.02 

3.61  ±0.01 

3.51  ±0.02 

3.21 

20  ss 

3.68  ±0.00 

3.59  ±0.03 

3.46  ±0.06 

3.29  ±0.01 

3.01  ±0.03 

ds 

3.92  ±0.04 

3.84  ±0.01 

3.72  ±0.00 

3.56  ±0.01 

3.26  ±  0.03 

(a) 


Buffer  Ionic  Strength  ImM] 


Figure  15:  Electrophoretic  mobility  of  10,  20  and  50  ss-  and  dsDNA  measured  in:  a)  1  pm  and  b)  100  nm 
deep  channels.  Data  points  corresponding  to  the  same  ionic  strength  have  been  plotted  slightly  apart  (in 
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the  x  direction)  for  the  ease  of  the  reader.  Error  bars  correspond  to  the  standard  deviation  of  all  trials  for  a 
given  experiment. 

In  all  cases,  as  the  ionic  strength  increases,  the  absolute  electrophoretic  mobilities  decrease, 
which  we  attribute  to  the  charged  solvent  effect  and  asymmetry  field  proposed  by  Manning,  and 
experimentally  confirmed  for  short  DNA  oligomers  (effect  (ii))  [17],  [26],  [34],  In  all  cases, 
ssDNA  absolute  electrophoretic  mobilities  are  smaller  than  dsDNA,  presumably  because  dsDNA 
has  a  higher  charge  density  [18].  It  is  also  interesting  to  note  that  recent  work  by  Stellwagen  and 
Stellwagen,  using  capillary  electrophoresis,  showed  that  the  separation  resolution  between 
ssDNA  and  dsDNA  was  inadequate  in  the  range  of  10-20  base  lengths,  and  that  10  base  ssDNA 
actually  had  a  higher  electrophoretic  mobility  than  10  bp  dsDNA  [35],  Their  experiments  were 
performed  in  0.5-2kV/m  applied  electric  fields,  with  a  detection  point  29.8  cm  downstream. 
Conversely,  our  experiments  were  performed  in  much  shallower  channels  (<  1pm),  which 
allowed  us  to  apply  ~  8x  larger  electric  field  strengths  without  significant  Joule  heating,  and  thus 
detect  at  shorter  distances  (1  cm).  This  fact,  combined  with  shorter  injection  plug  lengths,  allows 
us  to  minimize  diffusion  effects  and  greatly  increase  separation  resolution. 

To  further  understand  the  mechanisms  involved  in  DNA  separation,  we  investigate  the 
electrophoretic  mobility  difference  between  ss-  and  dsDNA  versus  buffer  ionic  strength  and 
strand  length,  shown  in  Figure  16,  in  both  1  pm  and  100  nm  channels.  The  dominant  parameter 
in  determining  the  difference  in  electrophoretic  mobility  between  ss-  and  dsDNA  is  the  shape  of 
DNA,  which  is  affected  by  salt-dependent  EDL  changes  around  the  DNA,  and  results  in 
differences  in  drag  (effect  (iii)).  The  DNA  shape  is  related  to  its  persistence  length,  Lp,  which  for 
dsDNA  (-45-50  nm)  is  larger  than  the  contour  length  of  any  of  the  strands  used  in  this  study,  for 
all  ionic  strengths  considered  [36].  Therefore,  in  this  study  only  ssDNA  will  exhibit 
conformational  changes  that  may  affect  electrophoretic  mobility.  However,  for  the  10  and  20 
base  ssDNA,  the  changes  in  Lp  are  too  small  to  induce  any  appreciable  change  in  conformation: 
indeed,  we  see  that  in  Figure  16  the  difference  in  mobility  for  the  20  bp  strands  remains 
approximately  constant.  Instead,  we  believe  that  the  downward  trend  we  see  in  Figure  16  for  50 
base  strands,  is  most  likely  due  to  these  conformational  changes.  When  ionic  strength  ranges 
from  10  to  60  mM,  the  ratio  of  strand  length  to  persistence  length  of  the  50  base  ssDNA 
increases  from  8.2  to  9.9  (whereas  it  is  only  a  shift  from  3.3  to  4  for  20  base  ssDNA,  and  1.6  to  2 
for  10  base  ssDNA).  These  ratios  correspond  to  a  46%  and  43%  mean  end-to-end  difference  for 
the  10  mM  case  and  60  mM  case  respectively  [37],  [38].  Thus,  we  conclude  that  as  the  ionic 
strength  increases,  the  50  base  ssDNA  tends  to  assume  a  more  compact  shape,  which  induces  a 
significant  increase  in  drag  and  results  in  a  higher  electrophoretic  mobility,  and  consequently,  in 
the  downward  trend  observed  in  Figure  15. 

For  10  bp  DNA,  we  hypothesize  that  the  mobility  trend  is  due  to  reaction  kinetics  within  the 
dsDNA  peak.  In  general,  the  observed  electrophoretic  mobility  of  dsDNA  is  dependent  on  the 
equilibrium  ratio  of  ss-  to  dsDNA:  when  the  equilibrium  ratio  is  large,  the  dsDNA  peak  is 
composed  of  mostly  ssDNA  and,  thus,  migrates  with  an  electrophoretic  mobility  that  is  closer  to 
the  electrophoretic  mobility  of  ssDNA.  For  example,  in  10  mM  buffer  the  melting  temperature 
(^=32.10)  is  closest  to  room  temperature  and,  therefore,  has  the  highest  equilibrium  ratio  of  all 
conditions.  In  this  case,  the  observed  dsDNA  electrophoretic  mobility  is  closest  to  the  ssDNA 
value,  and  results  in  poor  separation  resolution  (Figure  14a,  lower  trace).  As  Tm  increases 
(which  happens  as  the  ionic  strength  increases),  the  equilibrium  ratio  decreases:  separation 
resolution  improves,  as  the  dsDNA  peak  migrates  further  away  from  the  ssDNA,  and  its 
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eletrophoretic  mobility  gets  closer  to  its  true  value  (i.e.,  the  value  when  dsDNA  does  not 
disassociate).  Instead,  for  20  and  50  base  DNA,  at  all  ionic  strengths,  Tm  is  sufficiently  high  to 
result  in  extremely  low  equilibrium  ratios,  which  thus  have  negligible  effects  on  transport.  Note 
that,  as  the  DNA  elutes  down  the  separation  channel,  the  equilibrium  ratios  constantly  change 
due  to  separation  dynamics,  and  therefore,  the  calculated  melting  temperatures  are  only 
approximate. 


(b) 


lOOnm 


Buffer  Ionic  Strength  [mM] 


Figure  16:  Difference  of  electrophoretic  mobility  between  ss-  and  dsDNA  for  10,  20,  and  50  bp  DNA  in 
both  1  pm  (a)  and  100  nm  (b)  channels.  Data  shows  a  downward  trend  for  the  10  and  50  bp  (due  to 
changes  in  melting  temperature  and  persistence  length  respectively),  and  no  change  for  the  20  bp  ease. 
Error  bars  correspond  to  the  standard  deviation  of  all  trials  for  a  given  experiment. 

In  Figure  17,  we  investigate  the  effect  of  hydrodynamic  confinement  (effect  (iv))  more  closely 
by  plotting  the  percent  reduction  in  electrophoretic  mobility  between  a  1  pm  and  100  nm  deep 
channels.  First,  we  notice  that  the  behavior  of  ss-  and  dsDNA  is  qualitatively  very  similar.  This 
suggests  that,  at  least  for  DNA  smaller  than  20  bp,  confinement  into  nanometer  sized  channels 
affects  ss-  and  dsDNA  mobility  in  the  same  way.  Furthermore,  in  all  cases,  the  mobility  at  higher 
salt  is  reduced  by  about  4%.  We  attribute  this  reduction  to  an  increase  in  drag  due  to 
hydrodynamic  confinement. 


To  support  this  hypothesis,  we  analytically  solve  for  the  drag  on  a  spherical  particle  flowing 
between  two  walls,  using  the  analytical  expression  below  [39]: 


where  a  is  the  sphere  radius,  h  is  the  channel  height,  z  is  the  distance  from  the  center  of  the 
sphere  to  the  wall,  and  D0  is  the  free  solution  diffusivity. 
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10  Base  SS  10  Base  DS  20  Base  SS  20  Base  DS 


Figure  17:  Percent  difference  in  electrophoretic  mobility  between  1  pm  and  100  nm  channels  for  10  and 
20  bp  ss-  and  dsDNA  for  different  salt  concentrations.  Dashed  horizontal  line  represents  the  expected, 
theoretically  derived,  reduction  in  electrophoretic  mobility  due  to  increases  in  hydrodynamic  interactions 
with  the  wall  only  in  lOOnm  channels.  The  larger  differences  in  mobility  at  low  ionic  strengths  are  due  to 
electric  double  layer  interactions  between  the  walls  and  the  DNA.  Error  bars  correspond  to  the  standard 
deviation  of  all  trials  for  a  given  experiment. 

For  ssDNA,  we  approximate  a  as  being  equal  to  the  DNA  persistence  length.  This  way,  we  find 
that  the  hydrodynamic  drag  increase  for  ssDNA,  going  from  lum  to  lOOnm  channel,  is  about 
4.8%  (ssDNA  persistence  length  is  2.02  nm  at  60mM  [39]),  which  is  in  good  agreement  with  our 
data.  For  dsDNA,  the  spherical  approximation  is  geometrically  inaccurate,  as  these  strands  will 
tend  to  retain  a  rod-like  shape  at  all  salt  concentrations  studied  here.  At  the  same  time,  the  rigid 
cylinder  approximation  fails  to  give  meaningful  estimates  of  the  drag  coefficient,  because  it 
assumes  a  larger  aspect  ratio  than  DNA’s. 

Figure  17  also  shows  a  general  decreasing  trend  in  the  electrophoretic  mobility  difference 
between  1  pm  and  100  nm  channels  as  the  buffer  ionic  strength  increases,  coming  to  a  plateau  at 
larger  ionic  strengths.  As  ionic  strength  increases,  the  characteristic  size  of  the  EDL  or  the  Debye 
length,  Xd  ,  decreases  from  1.70  nm  at  lOmM,  to  1.03  nm  at  60mM,  and  any  electrostatic 
interactions  between  the  walls  and  DNA  become  negligible,  regardless  of  channel  height.  At  the 
lowest  ionic  strength,  the  EDL  accounts  for  3.4%  of  the  channel  height  for  the  lOOnm  deep 
channel  compared  to  0.34%  of  the  1pm  deep  channel.  We,  therefore,  hypothesize  that  the 
increase  in  electrophoretic  mobility  difference  between  DNA  in  different  channel  heights  at  low 
salt  concentration  is  due  to  electrostatic  interactions  between  the  DNA  and  the  EDL  (effect  (v)) 
that  cause  particles  to  move  towards  the  center  of  the  channel,  where  the  fluid  velocity  is  greatest 
[25],  This  reduces  the  time  to  arrival  at  the  detector  (because  electroosmotic  flow  is  larger  than 
electrophoretic  in  our  experiments),  making  the  calculated  electrophoretic  mobility  smaller  (see 
Equation  1).  Based  on  [25],  the  differences  in  overall  velocity  attributed  to  the  EDL  interactions 
for  a  particle  with  order  - 1  charge  would  be  on  the  order  of  one  percentage,  and  would  quickly 
become  negligible  as  Xo/h  decreases.  Thus,  we  believe  that  the  predominant  cause  for  the 
difference  between  the  mobilities  of  DNA  in  the  lum  and  lOOnm  channels  is  likely  due  to  such 
drag  effects.  A  full  theoretical/numerical  development  of  the  EDL  effects  on  both  the  DNA  and 
the  wall  is  beyond  the  scope  of  this  paper,  which  primarily  focuses  on  experimental 
measurements  with  heuristic  arguments  meant  to  support  plausible  hypothesis  explaining  our 
findings.  In  future  work,  we  will  pursue  more  detailed  theoretical  modeling,  which  so  far  we 
have  developed  only  for  the  simplified  case  of  spherical  particles. 
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A.2  Conclusions: 

Our  results  show  that  identification  of  short  (10,  20,  and  50  bases  long)  DNA  strands  using 
nanofluidic  channels  is  possible.  Although  the  initial  work  plan  included  the  investigation  of 
DNA  separation  and  identification  under  different  sample  conditions  (non  or  partially 
complementary  strands),  early  on  our  data  highlighted  the  complexity  of  DNA  separations  in 
channels,  and  the  various  effects  that  can  change  the  observed  behavior.  Therefore,  in 
consultation  with  Dr.  King,  it  was  decided  that  investigating  the  fundamental  science  underlying 
the  dynamics  of  DNA  in  micro  and  nanoconfmements  should  be  pursued  first,  as  a  prerequisite 
to  engineering  optimal  separation  platforms. 

What  we  found  is  that  buffer  ionic  strength  is  the  largest  factor  affecting  the  electrophoretic 
mobility  of  DNA,  but  we  also  showed  that  DNA  length,  conformation,  and  confinement  are 
important  variables.  Specifically,  we  showed  that  for  short  (lObp)  ss-  and  dsDNA,  one  needs  to 
consider  the  impact  of  melting  temperature  on  strand  mobility  and  separation  resolution;  while 
for  longer  strands  (50bp),  conformation  effects  play  an  important  role  in  determining  mobility. 
Finally,  we  showed  the  effects  of  confinement  on  the  absolute  electrophoretic  mobility  of  10  and 
20bp  ss-  and  dsDNA,  and  are  able  to  isolate  the  effects  of  hydrodynamic  confinement  and  EDL 
interactions.  These  results  are  the  basis  for  further  investigations  to  increase  the  sensitivity  and 
resolution  of  DNA  separations,  as  well  as  for  understanding  the  kinetics  of  DNA  with  low 
melting  temperatures.  Future  steps  will  include  further  improving  a  COMSOL  model  that  we 
started  to  develop  to  compare  experimental  results  with  theory,  and  that  will  prove  especially 
useful  in  understanding  the  coupling  between  separation  efficiency  and  DNA  hybridization 
dynamics. 
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